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ABSTRACT 


A  correlation  has  been  developed  between  the  inertial 
resistance  coefficient  and  reservoir  rock  properties  such  as 
porosity,  permeability  and  formation  resistivity  factor. 

Unlike  other  correlations,  this  is  dimensionally  consistent 
and  does  not  assume  the  Kozeny  model  for  porous  media. 

Experimental  determinations  on  twenty  reservoir  rock 
samples  have  shown  interrelationships  among  the  inertial 
resistance  coefficient,  porosity,  permeability,  formation 
resistivity  factor  and  tortuosity.  Data  on  porosity  and 
formation  resistivity  factor  indicates  that  their  product 
equals  approximately  four.  This  is  in  agreement  with  other 
published  data. 

A  method  has  been  devised  to  ensure  passage  of 
electric  current  through  brine  saturated  rock  samples  with 
fractured  ends,  without  developing  contact  resistances  at 
the  current  electrodes. 

Permeability  data  have  indicated  that  the  Klinken- 
berg  permeability  values  may  be  100%  to  200%  in  error  if 
proper  care  is  not  exercised  to  limit  the  flow  rate  to  the 
viscous  flow  region. 
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INTRODUCTION 


Fluid  flow  through  porous  media  may  be  described  by 
Darcy's  law  provided  the  flow  is  in  the  viscous  region  (1)* 

In  the  case  of  horizontal  laminar  flow,  the  flow  rate  varies 
directly  as  the  pressure  drop,  up  to  a  certain  point.  Beyond 
this  point,  both  the  proportionality  and  Darcy's  law  cease  to 
hold.  This  point  is  said  to  mark  the  boundary  between  the 
viscous  and  the  visco- inertial  regions  of  flow  (2)  . 

Horizontal  fluid  flow  in  both  the  viscous  and  visco- 
intertial  regions  may  be  described  by  the  equation: 

dP  9 

-  —  =  a  M  +  p  p  7^  (1) 

dX 

which  is  similar  to  an  equation  first  proposed  by  Forch- 
heimer(3).  In  order  that  this  equation.be  used  to  predict 
pressure  distribution  for  a  given  flow  rate,  it  is  necessary 
that  the  coefficients  a,  \±,  (3  and  p  be  known. 

The  coefficients  \x  and  p  represent  the  viscosity 
and  density,  both  well  known  fluid  properties.  As  a  result 
they  may  be  evaluated  either  by  measurement  or  estimation  by 
any  of  a  number  of  satisfactory  methods. 

The  coefficients  a  and  |3  are  referred  to  as  the 
viscous  and  the  intertial  resistance  coefficients.  They 
represent  rock  properties  which  are  dependent  upon  the 
geometry  of  the  porous  medium.  At  present  the  methods  avail¬ 
able  for  measuring  and  estimating  a  and  |3  are  not  entirely 
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satisfactory.  This  is  due,  in  part,  to  the  fact  that  a  means 
of  adequately  describing  the  geometry  of  a  porous  system  is 
not  available  and  that  the  exact  relationship  between  the 
geometry  and  the  coefficients  a  and  f3  is  not  yet  known. 

Since  a  is  the  inverse  of  the  rock  property,  per¬ 
meability,  it  may  be  evaluated  from  flow  test  using  the  Darcy's 
law,  provided  flow  is  restricted  to  the  viscous  region.  The 
inertial  resistance  coefficient  can  only  be  determined  experi¬ 
mentally  by  observing  flow  in  the  visco-intertial  region,  and 
employing  equation  (1),  or  a  modification  thereof.  The 
modified  form  of  this  equation  that  is  commonly  used,  may  be 
written  in  the  form  proposed  by  Green  and  Duwez(4). 


A  (Pz)  M  g  G 

— - =  a  +  (3  (-) 

2  ( ZRT)  LflG  M' 


(2) 


Although  this  method  is  adequate  for  research  work, 
it  is  considered  too  difficult  for  routine  core  analysis. 

To  circumvent  difficulty  various  attempts  have  been  made  to 
relate  (3  to  other  rock  properties  which  lend  themselves  more 
readily  to  evaluation  through  routine  core  analysis,  e.g. 
porosity,  permeability  etc.  The  first  such  attempt  was 
reported  by  Cornell  and  Katz (5)  who  correlated  (3  with  per¬ 
meability.  Janicek  and  Katz (6)  improved  the  correlation  by 
including  porosity  as  a  parameter.  Finally  Hamilton(7) 
extended  the  correlation  and  presented  a  new  expression  for 
(3  as  a  function  of  permeability  and  porosity. 
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An  examination  of  the  three  correlations  indicates 
that  they  are  dimensionally  inconsistent  and  that  they  are 
based  on  the  Kozeny  mod.el(8)  .  Since  a  dimensionally  incon¬ 
sistent  correlation  does  not  seem  reasonable  and  the  Kozeny 
model  does  not  admit  to  inertial  effects,  the  object  of  this 
study  was  to  examine  the  interrelationships  of  the  inertial 
resistance  coefficient,  permeability,  tortuosity,  porosity  and 
the  formation  resistivity  factor,  v/ith  the  view  of  developing 
a  correlation  free  of  the  two  above-mentioned  objections. 

This  study  was  restricted  to  an  examination  of  the 
data  of  previous  workers,  the  undamaged  rock  specimens 
studied  by  Hamilton(7)  and  a  group  of  ten  new  specimens. 

The  new  specimens  were  selected  to  provide  a  wide  range  of 
porosity  and  permeability. 
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THEORY 


Flow  Through  Porous  Media 

Flow  of  liquids  through  porous  media  was  first 
described  by  Darcy (9)  in  1856  by  the  equation 


q  =  -c 


dh 

dl 


(3) 


Although  this  equation  had  been  modified  and  used 
in  the  petroleum  industry  during  the  early  1930 's  it  was  first 
rigorously  examined  by  Hubbert(lO)  in  1940.  For  the  case  of 
horizontal  flow,  he  showed  that 


q 


K  dP 
|JL  dl 


(4) 


This  form  of  the  equation  was  found  to  be  satisfactory 
for  low  flow  rates  of  liquids.  However  in  order  that  the 
equation  be  applied,  K  and  \i  had  to  be  known.  For  application 
to  oil  reservoirs,  it  became  the  practice  to  subject  small 
rock  samples  to  flow  tests  and  apply  Darcy's  law  to  evaluate 
the  permeability,  K. 

Since  flow  tests  using  liquids  as  the  flowing  fluid 
proved  difficult,  it  seemed  appropriate  to  use  a  gas  instead 
of  a  liquid.  This  seemed  proper  since  it  had  been  shown  that 
K  was  a  property  of  the  porous  medium.  For  this  purpose  it  was 

shown  that  Darcy's  law  could  be  written  as 

K  d(P2 ) 

2|jlP„  dl 


q 


m 


(5) 


m 
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The  actual  application  of  this  method  for  obtaining 
K  was  found  difficult  since  even  for  low  flow  rates  K  varied 
with  pressure.  Klinkenberg( 11 )  examined  this  phenomenon  and 
showed  that,  whereas  with  liquids  there  was  no  appreciable 
slip  at  the  wall,  this  was  not  true  for  gases.  He  observed  that 
when  the  mean  free  path  of  the  molecules  of  a  gas  approaches 
the  size  of  the  flow  channel,  slip  at  the  wall  becomes  appreci¬ 
able  .  He  went  on  to  show  that  apparent  or  observed  permeability, 
Ka  varied  linearly  with  the  inverse  mean  pressure,  Pm. 

He  proposed  that  since, 


K 


a 


I<  ( 1  + 


b 


(6) 


m 


a  plot  of  Ka  versus  1/Pm  v/ould  yield  the  equivalent  liquid 

permeability  when  extrapolated  to  the  point  —  =  0 

Pm 

Subsequent  to  Klinkenberg ' s  observations,  it 
became  standard  practice  to  obtain  K  from  such  a  plot.  The 
permeability  value  so  obtained  w as  denoted  by  Kqq  and 
referred  to  as  the  Klinkenberg  permeability. 

The  procedure  consisted  of  conducting  flow  tests  by 
varying  the  upstream  pressure  and  maintaining  the  downstream 
pressure  at  atmospheric  conditions.  This  procedure,  often 
times  resulted  in  data  which  did  not  yield  a  linear  relationship 


Several  attempts  were  made  to  eliminate  the  slip 
effect  by  maintaining  it  at  some  fixed  value(12).  However, 
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Hamilton (7)  observed  that  although  such  a  practice  would  allow 
reproducibility,  different  values  of  mean  pressure  yielded 
different  permeabilities. 

Much  earlier,  it  had  been  observed  with  liquid  flow 

that  Darcy's  law  failed  at  high  flow  rates.  This  deviation 
from  Darcy's  law  was  attributed  to  turbulence  by  most  observers; 
whereas,  Hubbert(lO)  attributed  it  to  inertial  effects  and 
Houpeurt(13)  proposed  expansion  and  contraction  of  the  fluid 
during  its  flow  along  the  pore  channels,  as  the  cause. 

In  order  to  describe  flow  over  a  wide  range  of  flow 
rates  Green  and  Duwez(4)  suggested  the  equation 

-  —  =  a^xv  +  (3pv2  (7) 

dL 


which  is  similar  to  an  equation  first  proposed  by  Forchheimer ( 3) . 
For  the  particular  case  of  gas  flow,  they  reasoned  that  if 
pv  -  G  »  constant,  then  the  equation  could  be  written  as 


or  as 


M  gc  d(P2 )  2 

- -  _  a  li,  G  +  (3  G  (8) 

2 ( ZRT)  dL 


A  (P2)  M  g  G 

- c_  =  a  +  (3  (-) 

2L  (ZRT)  [l  G 


(9) 


in  the  integrated  form. 

A  similar  result  was  obtained  by  Cornell  and  Katz (5). 
They  considered  the  deviations  from  Darcy’s  law  to  be  the 
result  of  turbulence.  Using  a  Kozeny  type  equation(14)  in 


- 


-  7 


terms  of  effective  pore  diameter  (D  )  and  tortuosity  r,  they 

.L j 

developed  the  equation 


L 


32F 


2ZRT 


KlV  M 


G 


32F1,b  2ZRT 

*  0.5* 

K.Dn  M 

Z  xh 


G 


(10) 


The  factors  and  K2  were  termed  as  geometrical 
constants,  and  had  to  be  determined  for  each  porous  medium. 

The  constant  ly  v/ac  closely  related  to  the  Kozeny  constant ( 14) 
and  was  assumed  to  equal  0.5.  The  factor  I<2  was  correlated 
with  the  relative  roughness  £  ,  as  in  the  case  of  pipes.  They 
assumed  that  some  sort  of  correlation  must  exist  between  I<2 
and  £  ,  since  the  latter  is  characteristic  of  a  porous  medium. 

A  comparison  of  equations  (9)  and  (10)  shows  them 
to  be  identical,  provided  that, 


a 


32F 

< - J> 

kide 


32F 


1.5 


and 


(3 


K  D  (fi 
2  E 


0.5 


(11) 


Initially  these  equations  were  of  academic  interest 
only.  However  an  examination  of  certain  natural  gas  reser¬ 
voirs  showed  that  at  high  delivery  rates,  flow  in  the  reservoir 
was  in  the  visco-inertial  region.  These  equations  then  became 
important  in  the  study  of  such  reservoirs.  In  order  to  make 
such  studies  values  of  the  inertial  resistance  coefficient 
p>  had  to  be  known.  Although  (3  could  be  obtained  experimentally, 
techniques  were  sought,  for  estimating  p>  from  other  parameters. 
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In  the  earliest  such  attempt  Cornell  and  Katz (5) 

proposed  the  correlation 


P 


4.11  x  1010 
j^l  •  3  3 


(12) 


Subsequently  Janicek  and  Katz (6)  re-examined  the 

data  of  Cornell  and  observed  that 

F0  =  4  (13) 

9.411  x  1013 

By  assuming  -  0.5,  and  noting  that  F0  -  4  and  K  ~ - — — - - 

a 


they  showo  d  th  a  t 


8  ( FK) °  *  b 


D 


E 


9.7  x  10 


6 


(14) 


by  virtue  of  equation  (11). 

Sensing  that  the  factor  K2  was  related  to  0  and  I< 
they  plotted  K2  versus  K/0  and  expressed  the  best  fit  by 


-2  ,K  0,75 

K0  =  2.82  x  10  (~) 

2  0 


(15) 


Finally  substituting  for  Dp,  F  and  K2  in  equation  (11)  they 

obtained  the  expression 


9 

5.5  x  10 

P  =  - - —  (16) 

k1.25  0O.75 


which  they  proposed  as  a  general  correlation ( 15) . 

Hamilton (7)  observed  that  the  Janicek  and  Katz 
correlation  did  not  agree  well  with  experimental  data.  He 
therefore  obtained  additional  data  and  on  the  basis  of  these 
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along  with  those  of  Cornell  he  proposed  that 


K. 


F 


( 17 ) 


which  when  substituted  into  equation  (11)  yielded  the  expression 


6.6  x  10' 

K1 .13  ^2 .76 


(18) 


He  proposed  this  as  a  general  correlation. 


E 1 ec t r i c a 1  P r  op e r t i o s  o f  P or ou s  Mo d  1  a 

The  electrical  properties  of  a  rock  depend  upon  the 
geometry  of  the  voids  and  the  fluids  which  fill  them.  The  use 
of  such  properties  to  characterize  reservoir  rocks  was  first 
suggested  by  Archie (16).  Amongst  other  things  he  defined  the 
formation  resistivity  factor  as 

Resistivity  of  a  saturated  medium 
Resistivity  of  the  saturating  solution 

This  factor,  is  a  parameter  which  reflects  the 
influence  of  pore  geometry  on  electrolytic  conduction  through 
the  medium.  It  in  effect  combines  porosity  and  tortuosity 
into  one  parameter. 

Since  resistance  of  a  saturated  medium  may  be 
expressed  as 


(0A) 


Pb 


(19) 
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and.  since  the  resistance  of  a  congruent  volume  of  the  saturating 
fluid,  may  be  written  as 


R 


L 


(20) 


b 


then  by  definition 


F 


(Pb  Lc)/(0A) 


(pb  L)/A 


(21) 


However,  the  square  of  the  ratio  of  the  tortuous  average 
length  of  the  pore  channels  traversed  by  the  current  flov/ing 
from  one  to  the  other  of  the  parallel  planes  within  the 
saturated  medium  to  the  distance  of  the  two  planes,  has 
been  defined (17)  as  the  tortuosity  of  the  medium.  Therefore, 
the  formation  resistivity  factor  may  be  expressed  as 


0.5 


T 


F 


Wyllie  and  Rose (18)  concluded  that  the  hydraulic 


and  electrical  tortuosities  in  a  porous  medium  are  extra¬ 
ordinarily  close  and,  for  all  purposes  could  bo  considered 
identical . 

Archie (16)  correlated  the  observed  formation 
resistivity  factor  with  porosity  and  permeability.  He 


' 


-  11 


.suggested  that  the  correlation  with  </  was 
and  that  the  formation  resistivity  factor 
as 


the  better  correlation 
could  be  expressed 


F 


(23) 


where  m  is  the  cementation  factor.  Ha  further  reported  that 
m  probably  ranged  from  1.8  to  2.0  for  consolidated  sand¬ 
stones  . 

Winsauer(19)  reported  a  similar  relationship  based 
on  correlations  of  their  data  from  a  large  number  of  sand¬ 
stone  cores.  This  equation,  commonly  referred  to  as  the 
Humble  equation (20)  is 

F  =  0.62  0~2*lj  (24) 


At  present  it  is  best  to  say  that  the  formation 
resistivity  factor  is  a  function  of  porosity  and  some  unknown 
effect  of  the  complex  internal  geometry  of  the  porous  medium. 
Win saucr  and  Shear in  (19)  and  Wyllie  and  Spangler (17)  have 
summarized  the  pertinent  literature  on  this  aspect. 


12 


EXPERIMENTAL  EQUIPMENT 

It  was  decided  to  obtain  data  for  the  inertial 
resistance  coefficient,  porosity,  permeability  and  formation 
resistivity  factor  on  ten  cores  used  by  Hamilton (7)  as  well 
as  an  additional  set  of  ten  selected  core  samples.  To 
accomplish  this,  various  pieces  of  equipment  were  necessary. 

Sample  Preparation  Equipment 

The  equipment  used  to  prepare  the  rock  specimens 
for  porosity  and  flow  measurements  was  identical  to  that 
used  by  Hamilton(7) .  In  order  to  prepare  the  specimens  for 
resistance  measurements,  a  gram-atic  balance,  a  vacuum  jar, 
a  mechanical  vacuum  pump  and  a  de-aerated  50,000  p.p.m.  brine 
solution  were  used* 

Porosity  Measurement  Apparatus 

A  modified  Oilwell  Research  Porosimeter ( 21)  using 
helium  was  employed  to  obtain  data  on  grain  volume.  An 
electric  Mercury  Pycnometer ( 18)  was  used  to  obtain  data  on  the 
bulk  volume  of  the  samples.  Both  these  pieces  of  equipment 
were  calibrated  prior  to  use  by  means  of  six  solid  blanks 
of  known  dimensions. 
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Permeability  Measurement  Apparatus 

Several  method  for  holding  small  rock  samples  for 
the  purpose  of  measuring  permeability  are  described  in  the 
literature (5  )(  22)(  23).  Hamilton(7)  modified  the  method 
described  by  Johnston  and  Taliaf erro( 2 3)  and  employed  a  core 
holder  machined  from  mild  steel  stock.  In  this  work,  the  same 
core  holder  was  employed.  A  line  diagram  of  the  permeability 
measurement  apparatus  is  shown  in  figure  1. 

The  core  holder  consisted  of  two  pressure  taps  in 
the  body  (upstream  and  downstream),  and  two  end  plates  with 
0  ring  seals.  Each  end  plate  was  provided  with  a  threaded 
pressure  fitting,  and  an  Iron-Constantan  thermocouple. 

A  sleeve  was  made  for  each  core  using  an  extra  long 
No.  26  rubber  stopper.  Each  stopper  was  drilled  to  a  diameter 
0.1  inch  less  than  the  diameter  of  the  core  for  which  it  was 
intended.  This  provided  for  a  squeeze  fit  of  the  core  into 
the  sleeve.  The  sleeve  containing  the  core  was  placed  into 
the  tapered  section  of  the  core  holder,  followed  by  a  sliding 
metal  follower  and  a  SKF  No.  51105  thrust  bearing.  These 
were  held  in  place  by  screwing  down  a  metal  follower. 

Nitrogen  was  used  as  the  flowing  fluid.  The  inlet 
flowing  pressure  was  regulated  at  the  gas  cylinder  by  means 
of  a  3000/200  p.s.i.  regulator  for  low  pressures  and  a 
3000/3000  p.s.i,  regulator  for  pressures  in  excess  of  200  p.s.i. 
From  the  regulator  the  gas  was  passed  through  a  0.25  inch 


0-2000  psi  t  0-250  psi 

0-  100  psi 


14 


l 


15 


diameter  soft  copper  tubing  and  then  admitted  to  the  core 
holder  through  a  0,125  inch  stainless  steel  tubing.  The 
copper  tubing  was  immersed  in  a  glycol  bath,  maintained  at  a 
constant  temperature  by  means  of  a  heating  and  refrigeration 
unit,  so  as  to  control  the  temperature  of  the  flowing  fluid. 
The  sensitivity  of  the  bath  was  -0.5°F. 

Two  stainless  steel  needle  valves  were  arranged  in 
parallel  on  the  downstream  side  of  the  core  holder  to 
provide  for  the  application  of  back  pressure,  if  and  when 
desired. 

Pressure  readings  were  obtained  by  means  of  water 
and  mercury  manometers  for  pressures  up  to  30  p.s.i,,  a 
Seeger‘s  gauge  for  the  range  30  -  100  p.s.i.,  a  250  p.s.i. 
Heise  gauge  for  the  range  100  -  250  p.s.i.  and  2000  p.s.i. 
Heise  gauge  for  the  range  250  -  2000  p.s.i.  All  these  in¬ 
struments  were  calibrated  prior  to  use. 

Flowing  gas  temperatures  at  the  upstream  and 
downstream  ends  were  measured  by  means  of  the  previously 
mentioned  Iron-Constant an  thermocouples,  and  a  Leeds  and 
Northrup  potentiometer . 

Gas  flow  measurements  were  made  using  an  electric 
stop  watch  and  various  volume  measuring  instruments,  depending 
upon  the  rate  of  flow,  as  follows: 

a.  Selas  Porcelain  Microfilter  tube  0<Q^<0.5 

0.5$Q<2 


b.  Gas  Buret 
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c,  Wet  Test  Meter  2<Cq^20 

d.  Diaphragm  Meter  20  ^120 

where  Q  is  the  rate  of  flow  in  ft3/hr. 

Resistance  Measurement  Apparatus 

An  apparatus  for  measuring  resistance  of  small  rock 
samples  is  described  in  the  literature  by  Rust(24).  A 
modified  version  of  his  apparatus  was  found  more  appropriate 
for  this  investigation.  The  apparatus  consisted  of  a  lucite 
housing,  split  into  similar  halves,  which  could  accomodate 
the  rock  sample.  Two  sets  of  six,  solid  silver,  platinum 
tipped,  potential  electrodes,  under  spring  tension  were 
seated  in  the  housing  in  such  a  way  that  when  the  core  was 
in  place,  the  electrodes  were  distributed  along  the  circum¬ 
ference  of  the  core  sample.  Two  brass  discs  were  used  as 
current  electrodes.  They  were  arranged  to  make  contact  with 
the  two  ends  of  the  core  sample  through  brine  saturated  felt 
pads.  The  pads  were  threaded  with  thin  silver  wire  to  ensure 
electrical  contact  between  the  brass  discs  and  the  core. 

The  lucite  housing  and  the  brine  saturated  core 
were  housed  in  a  lucite  box,  partially  filled  with  50,000 
p.p.m.  brine  to  minimize  evaporation  while  the  test  was  in 


progress , 


Figurp  2.  Resistance  Measurement  Cell 
(after  Rust(24)) 
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The  wiring  of  the  apparatus  was  of  the  permanent 

type  as  shown  by  the  line  diagram  in  Figure  2a.  Two  selector 
switches  provided  for  resistance  measurements  by  either  the 
two  electrode  or  the  four  electrode  method. 

Dunlapfs(25)  check,  for  contact  resistance  developed 
at  the  current  electrodes,  was  incorporated  by  including  a 
series  resistor  equal  in  value  to  the  injjut  impedance  of 
the  vacuum  tube  voltmeter.  The  voltmeter  used  was  a 
Hewlett  Packard  model  400  AB.  The  input  line  was  shielded 
so  as  to  minimize  the  effect  of  background  noise. 
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Variable 

A.C 

Source 


50,000  U  ,  adjustable  resistor 
PR  100  &  ;  precision  resistor 

VTV  A-C  vacuum  tube  voltmeter 

ccR  contact  check  resistor;  10  megi^ 
PP  Potential  probes 

S^  Selector  switch  1 

S2  Selector  switch  2 


Circuit  Diagram  for  Resistance  Measurements 
(after  Rust(24)) 


Figure  2-a 
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EXPERIMENTAL  PROCEDURE 


In  addition  to  ten  cores  used  by  Hamilton (7 ) ,  ten 
new  cores  were  selected  for  observation.  The  selection  was 
such  as  to  provide  a  wide  range  of  porosity  (2%  -  27,5%)  and 
permeability  (0.000125  md  to  1160  mds) . 

Sample  Preparations 

Each  small  cylindrical  reservoir  rock  sample  was 
cut  parallel  to  its  bedding  plane  using  a  one-inch  diamond 
drill,  and  was  limited  in  length  to  three  inches.  Since  it 
was  observed  that  the  end  pores  became  plugged  due  to  the 
sawing  action,  new  end  faces  were  produced  by  fracturing,  using 
the  device  described  by  Hamilton(7),  The  fractured  cores 
were  approximately  1.5  inches  in  length  and  one  inch  in 
diameter . 

The  physical  dimensions  of  the  cores  were  obtained 
by  means  of  a  micrometer  caliper.  The  effective  diameter  of 
each  core  was  taken  to  be  the  average  of  five  readings  and  the 
effective  length  as  the  average  of  seven  readings. 

The  cores  were  extracted  by  means  of  a  Soxhlet 
apparatus  using  benzene  as  the  solvent.  After  an  extraction 
period  of  about  18  hours,  the  cores  were  allowed  to  soak 
over  night  in  fresh  benzene.  If,  after  this  period,  the 
benzene  was  found  to  be  clear,  it  was  inferred  that  the 
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extraction  was  complete.  If,  however,  the  benzene  was  oil- 
stained,  the  procedure  was  repeated.  Subsequent  to  satis¬ 
factory  extraction,  the  cores  were  dried  in  a  constant  tempera¬ 
ture  oven  at  200 °F  for  about  12  hours.  The  cores  were  allowed 
to  cool  to  room  temperature  and  were  then  stored  in  a  dessicator 
until  used. 

Porosity  Measurements 

The  data  on  matrix  volume  of  the  cores  were  obtained 
using  the  modified  Oil  Well  Research  Porosimeter ( 21) .  The 
data  on  bulk  volume  were  obtained  using  an  electric  mercury 
Pycnometer ( 18) .  Six  steel  blanks  and  six  lucite  blanks  of 
known  dimensions  were  used  to  calibrate  the  porosimeter  and 
the  pycnometer,  respectively.  From  these  data,  the  porosity 
values  were  computed  and  recorded  in  Table  1.  The  observations 
were  made  according  to  standard  procedures,  therefore,  no 
further  details  are  given  here. 

Permeability  Measurements 

Each  core  housed  in  its  rubber  sleeve  was  placed 
in  the  tapered  section  of  the  core  holder.  The  sleeve  was 
slightly  greased  on  the  outside  to  facilitate  loading  and  to 
prevent  fluid  bypass  during  test.  The  follower  was  screwed 
down  tightly  and  the  end  plates  were  assembled  by  tightening 


22 


Table  1 


Physical  Properties 

i  and  Porosity 

Data 

Source 

Core 

No. 

Formation 

Rock  Type 

Core 

Diameter 
( inch) 

Effective 
Length 
( inch) 

Hamilton 

3 

Cardium 

Sandstone 

0.998 

1.490 

Hamilton 

4 

Cardium 

Sandstone 

0.997 

1.461 

Hamilton 

5 

Cardium 

Sandstone 

0.996 

1.462 

Hamilton 

7 

Sturgeon  Lake 

Limestone 

0.996 

1.498 

Hamilton 

8 

Berea 

Sandstone 

0.996 

1.474 

Hamilton 

10 

Blairmore 

Basal  Quartz 

1.015 

1.445 

Hamilton 

11 

Blairmore 

Basal  Quartz 

1.015 

1.341 

Hamilton 

12 

Viking 

Sandstone 

1.015 

1.414 

Hamilton 

13 

Viking 

Sandstone 

1.016 

1.467 

Hamilton 

14 

Belly  River 

Sandstone 

1.015 

1.480 

C  &  G  Lab** 

4X 

Cardium 

Sandstone 

0.9820 

1.4931 

C  &  G  Lab . 

20 

Belly  River 

Sandstone 

0.9860 

1.4760 

C  Sc  G  Lab, 

22 

* 

Sandstone 

0.9934 

1.5292 

C  Sc  G  L  ab . 

24 

* 

Sandstone 

0.9996 

1.5746 

C  Sc  G  Lab, 

31 

Belly  River 

Sandstone 

0.9880 

1.5100 

C  Sc  G  L  ab , 

38 

* 

Sandstone 

0.9934 

1.5770 

C  Sc  G  L  ab , 

45 

Taber 

Sandstone 

0.9900 

1.526 

C  Sc  G  Lab, 

46 

Taber 

Sandstone 

1.0000 

1.5601 

C  Sc  G  L  ab , 

53 

Taber 

Sandstone 

0.9934 

1.5598 

C  Sc  G  Lab, 

61 

Taber 

Sandstone 

0.9996 

1.5922 

*  Unknown 

**  The  Chemical 

and  Geological  Laboratories 

,  Edmonton 

Table  1  -  continued 


Core 

No. 

Bulk  Volume 
(cc) 

Pore  Volume 
(cc) 

0 

% 

Aiea  Effective  Length  (L) 

m . 

Area  (A) 

3 

19.094 

3.593 

18.82 

0.782 

1.9054 

4 

18.690 

3.945 

21.11 

0.781 

1.8707 

5 

18.664 

3.202 

17.16 

0.779 

1.8768 

7 

19.150 

2.605 

13.60 

0.780 

1.9205 

8 

18.806 

3.646 

19.39 

0.778 

1.8946 

10 

19.165 

3.82  5 

19.96 

0.810 

1 .7  840 

11 

17  .793 

3.628 

20.39 

0.810 

1.6556 

12 

18.7  57 

1.842 

9.82 

0.810 

1.7457 

13 

19.472 

2.002 

10.28 

0.810 

1.8111 

14 

19.629 

3.701 

18.85 

0.810 

1.8272 

4X 

18.82 

4.17 

22.18 

0.7574 

1.9712 

20 

18.28 

0.42 

2.30 

0.7636 

1.9329 

22 

18.95 

4.60 

24.33 

0.7775 

1.9668 

24 

19.77 

4.28 

20.91 

0.7848 

2.0064 

31 

18.50 

1.08 

5.84 

0. /667 

1.9695 

38 

19.32 

4.50 

23.36 

0.7775 

2.0026 

45 

18.10 

3.90 

21.55 

0.7698 

1.9823 

46 

19.68 

3.73 

18.10 

0.7854 

1.9864 

53 

19.47 

5.05 

22.06 

0.7775 

2 .0062 

61 

19.87 

5.69 

27.43 

0.7848 

2.0288 

24 


the  five  threaded  steel  bolts  provided  for  that  purpose. 

The  core  holder  assembly  and  the  connections  were 
tested  for  gas  leakage  by  pressuring  to  1800  p.s.i.  and  using 
standard  soap  solution.  The  soap  solution  was  examined  for 
bubbles  by  means  of  a  magnifying  lens.  Any  leaks  that  were 
detected,  were  removed  by  coating  the  area  with  Lepage's  glue. 

Nitrogen  at  a  regulated  pressure  and  temperature 
was  admitted  into  the  core.  The  downstream  end  pressure  and 
the  flow  rate  was  measured  using  instruments  previously 
mentioned.  In  accordance  with  the  Klinkenberg  theory (11)  ,  the 
flow  rate  was  limited  to  the  viscous  flow  region  by  applying 
the  desired  amount  of  back  pressure.  The  temperatures  of 
the  gas  at  the  upstream  and  the  downstream  ends  were  observed 
and  recorded. 

Several  readings  were  taken  at  various  mean  pressures, 
and  the  corresponding  flow  rates  were  noted.  Three  flow  rate 
readings  were  taken  for  each  data  point  to  make  sure  that 
stabilized  flow  was  attained.  Whenever  the  upstream  pressure 
was  changed,  a  stabilization  period  of  from  15  -  20  minutes 
was  provided. 

V i SCO- Inertia .1  Flow  Measureme nts 

Since  visco-intertial  flow  was  desired,  back  pres¬ 
sure  was  not  applied.  As  in  the  case  of  the  permeability 
measurements,  three  readings  were  taken  for  each  data  point. 
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Sample  Prep ar at ions  for  Resistance  Measurements 

The  clean  dry  cores  were  first  weighed.  Then  they 
were  evacuated  to  0.15  mm  of  mercury  by  means  of  a  vacuum 
pump.  A  de-aerated  50,000  p.p.m.  sodium  chloride  solution 
(prepared  by  adding  pure  sodium  chloride  to  freshly  boiled 
distilled  water)  was  slowly  admitted  into  the  vacuum  jar 
until  all  the  cores  were  submerged  in  the  solution.  The 
cores  were  then  allowed  to  remain  in  the  submerged  condition 
for  12  hours. 

The  saturated  cores  were  weighed.  This  was  done 
by  placing  the  saturated  cores  on  a  saturated  felt  pad  so 
as  to  minimize  the  loss  of  brine  by  drainage.  Using  these 
data,  the  saturations  were  calculated,  as  given  in  Table  4. 

Resistance  Measurements 

A  schematic  wiring  diagram  of  the  resistance 
measuring  circuit  is  given  in  figure  2a.  With  the  core 
sample  mounted  in  the  lucite  housing,  the  current  flowing 
through  the  core  was  set  at  one  milliampere  by  adjusting 
the  voltage  drop  across  the  wire  wound  precision  resistor. 
Manipulation  of  switches  I  and  II  permitted  voltage  drop 
measurements  to  be  made,  either  across  the  entire  core  sample 
(two  electrode  method)  or  across  a  portion  of  the  sample 
(the  four  electrode  method) * 
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The  current  electrodes  were  checked  for  contact 
resistance  by  Dunlap^s  method(25)*  If  and  when  such  resis¬ 
tance  was  found  to  be  present,  it  was  removed  by  re- saturating 
the  felt  pads  and  re-making  of  the  current  electrode  contacts. 
Subsequent  to  the  removal  of  contact  resistance,  three  readings 
were  taken  on  each  core  by  each  of  the  two  methods. 

The  resistivity  of  the  brine  was  measured  periodically 
by  means  of  a  standard  conductivity  bridge. 
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DISCUSSION  OF  EXPERIMENTAL  RESULTS 


The  data  on  porosity,  permeability,  formation  resis¬ 
tivity  factor  and  tortuosity  were  obtained  using  the  techniques 
outlined  elsev/here  in  this  work. 

The  permeability  data  were  obtained  by  applying  back 
pressure  to  limit  the  flow  rates  to  the  viscous  region.  This 
was  found  necessary  because  anomalous  values  of  Klinkenberg 
permeability,  which  were  at  least  100%  lower  than  the  true 
Klinkenberg  permeability,  were  obtained  if  the  required  amount 
of  back  pressure  was  not  applied. 

It  is  a  common  practice  to  plot  a  set  of  1/P 

values  against  corresponding  values  of  apparent  permeability, 

K  ,  and  to  extrapolate  the  straight  lino  portion  to  =  0, 

a  '  pm 

to  give  the  Klinkenberg  permeability.  It  is  also  usual  to 
regard  the  straight  line  portion  of  the  plot  to  indicate 
absence  of  the  slip  effect,  and  any  deviation  from  the 
straight  line  to  indicate  the  slip  phenomenon . 

This  study  has  indicated  that  this  criterion  does 
not  always  hold  ture .  Data  on  core  No.  53  have  been 
included  here  to  elaborate  this  statement.  A  comparison  of 
Figures  3  and  4  shows  that  the  Klinkenberg  permeability 
evaluated  from  data  obtained  without  application  of  back  pressure 
is  175  rads,  and  375  rads  when  evaluated  from  data  taken  when  back 
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Figure  3.  Klinkenberg  Permeability  Plot 
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Figure  4.  Klinkenberg  Permeability  Plot 
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pressure  was  applied.  In  Figure  3  a  straight  line  portion  is 
indicated  for  values  of  —  ^  0.04  (psia)  “  .  This,  however, 

does  not  necessarily  mean  that  the  data  in  this  range  are  in 
the  viscous  region.  In  fact  it  is  very  likely  that  these 
points  are  not  in  the  viscous  region  since  points  of  equal 
1/P  value  in  Figure  4  yield  higher  permeability  values. 

It  was  therefore  decided  to  use  only  viscous  flow 
readings  for  evaluating  the  Klinkenberg  permeability.  This 
was  achieved  by  taking  numerous  data  and  selecting  only 
those  which  were  in  the  viscous  region, 

2 

Katz  et  al(15)  suggested  that  a  plot  of  A(P  )  versus 
Q  should  be  adequate  to  represent  the  flow  data  in  both  the 
viscous  and  the  visco-inertial  regions.  The  straight  line 
portion  of  the  plot  represented  viscous  flow  and  any 
deviation  from  the  line  indicated  the  visco-inertial  flow 
region,  A  plot  of  this  type  was  made  for  each  core.  The 
viscous  flow  region  was  determined  in  each  plot  as  is  shown 
for  core  53  in  Figure  5. 

The  permeability  data  were  also  evaluated  using 
the  a  values  obtained  from  the  flow  data  in  the  visco-inertial 
region.  Figure  6  shows  the  permeability  value  so  obtained. 

A  comparison  of  the  values  obtained  by  these 
two  techniques  is  shown  in  Table  2. 
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Figure  5.  Flow  Regime  Plot 


32 


36  W  (  d)V 


GO¬ 

'S 

§ 

O 


Figure  6.  Flow  Data  Plot  to  Calculate 


Summary  of  Porosity,  Permeability  and  Inertial  Resistance  Coefficient  Data 
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Figures  1A  to  18A  in  Appendix  II  show  the  permeability 
values  obtained  by  the  above-mentioned  techniques. 

In  order  to  arrive  at  a  general  correlation  between 
the  inertial  resistance  coefficient  and  the  reservoir  rock 
properties,  Hamilton's  data(7)  and  Cornell's  data(14)  were 
also  included  in  this  work.  An  examination  of  Figures  12, 

13  and  14  shows  that  scatter  of  Cornell's  data  is  more  than 
either  Hamilton's  work  or  this  work.  There  is  reason  to 
question  the  accuracy  of  Cornell’s  data. 

In  order  to  obtain  values  of  porosity  Cornell 
saturated  his  samples  with  distilled  water.  Such  a 
technique  will  yield  results  of  questionable  accuracy  for  two 
reasons.  First,  porosity  determination  by  liquid  saturation 
is  inherently  inaccurate  since  complete  saturation  is  dif¬ 
ficult  to  attain  as  may  be  seen  from  Table  4,  Appendix  1. 

Second,  it  is  a  well  known  fact  that  the  introduction  of 
distilled  water  into  a  natural  porous  medium  will  cause 
swelling  of  argillaceous  impurities,  thereby  altering  both 
the  porosity  and  permeability. 

The  permeability  values  used  by  Cornell  were  cal¬ 
culated  from  a  values  only  and  were  not  checked  by  the 
Klinkenberg  method. 

An  examination  of  Equations  12,  16  and  18  representing 
the  correlations  proposed  by  Cornell  and  Katz(5),  Janicek  and 
Katz(6)  and  Hamilton(7)  showed  that  they  are  based  on  a  modifica¬ 
tion  of  Kozeny ' s  model (14)  for  porous  media.  The  Kozeny 
model  assumes  a  rectilinear  system  of  pore  channels,  and 
thus  cannot  account  for  the  effect  of  inertial  forces. 
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which  develop  because  of  a  curvilinear  system  of  pore  geometry 
(10) (13).  Also,  those  correlations  are  dimensionally  incon- 
sistent  because  p  (ft  )  cannot  be  equated  to  K  (ft  ),  In 
view  of  these  objections  the  very  basis  of  the  correlations 
was  re-examined. 

The  internal  resistance  coefficient  (3  is  known  to 
depend  on  the  pore  geometry  of  a  porous  medium.  Parameters 
such  as  0,  K,  F,  S^,  t ,  etc.  characterize  the  pore  geometry. 
Intuitively  therefore,  a  relationship  must  exist  between  (3 
and  the  reservoir  parameters.  This  relationship  may  be 
expressed  mathematically  as 

p  =  f^ 'p -.-Te  geometry) 

where  f ^  denotes  some  function  . 
and 

pore  geometry  =  (0,  K,  F,  t,  Sa,  etc) 

where  is  some  function. 

Therefore, 

P  =  f1  [f2  (0,  K,  F,  t,  Sft,  etc) 

=  f3  (0,  K,  F,  t,  SA,  etc)  (25) 

where  £■$  is  an  explicit  function, 
or 

f4  (P,  0/  K,  F,  t,  SA,  etc)  = 
where  f4  is  an  implicit  function. 


0 


(26) 
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Also,  it  has  been  shown  by  previous  workers ( 17 )( 19 )( 20) 

that 

r  ="  (F,  0)  where  f^  is  a  function  (27) 

and 

F  =  fg  (0)  v/here  is  a  function  (28) 

During  this  investigation,  no  attempt  was  made  to 
collect  data  on  internal  surface  area,  and  tortuosity  by  a 
means(19)  independent  of  the  formation  resistivity  factor. 

Therefore,  equation  26  was  reduced  to 

f?  O,  0,  K)  -  0  (29) 

where  f^  is  an  implicit  function  relating  the  variables  (3,  0 
and  K. 

0 . 5 

A  plot  of  log  (3  versus  log  (K)  *  indicated  relation¬ 
ship  as  shown  in  Figure  8.  However ,  a  plot  of  log  (3  versus 
log  0  also  showed  a  relationship,  as  indicated  in  Figure  9. 

In  order  that  the  proposed  correlation  be  dimen¬ 
sionally  consistent  it  was  observed  that  the  permeability 
2  0  5 

K(ft  )  should  appear  as  (K) in  a  grouping  of  0,  K,  f3,  and  F. 

A  factor  f  was  defined  such  that 
s 

(9.413  x  1013)°’b 

fs  =  — — — — - — (30) 

P-  (K)0*5 

v/here  the  numerator  is  a  conversion  factor  to  convert  per¬ 
meability  from  ft^  to  mds . 
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Log  /K  Figure  8 
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Log  ^  (  fraction  ) 

Figure  9 
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By  applying  the  method  of  averages(28)  to  Cornell's, 
Hamilton's  and  the  data  obtained  during  this  study,  a  plot  of 
log  f  versus  log  0  resulted  in  the  equation 

fg  =  60  05  * 07  (31) 

as  shown  in  Figure  11.  The  substitution  of  this  value  for 

f  into  equation  (30)  results  in  the  following  expression  for 
s 

the  inertial  resistance  coefficient. 


P  =  (1.617  x  105)  K-0*5  0~5'07  (32) 


This  expression  is  proposed  as  a  general  correlation. 
Although  it  can  be  used  in  this  form,  the  same  correlation  may 
be  represented  graphically  as  indicated  in  Figure  12.  In  the 
event  that  equation  (32)  is  used  to  calculate  (3,  the  rounding 
of  the  exponent  of  0,  from  5.07  to  5.0  seems  quite  permissable 
as  shown  in  Figure  13. 

Further  improvement  in  the  proposed  correlation  was 
attempted  by  incorporating  the  formation  resistivity  factor 
in  the  grouping  {0  K  *  ) .  Several  groupings  with  different 


indices  of  F  were  tried.  The  grouping 


Q*5-  (KF)°-5J 


appeared 


to  fit  the  data  best.  A  plot  of  log  (3  versus  log 


[_05.  (KF)°-5J 


is  shown  in  Figure  14.  It  is  important 

to  mention  here  that  this  correlation  shows  less  scatter 

♦ 

i 

than  any  other  correlation  presented  so  far. 
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Log  <f>  (fraction) 


Figure  11 


-  i  '  - 


0  (ft'1) 


'I 

«*  .u 


K  (mds) 


Figure  12.  Correlation  of  Intertial  Resistance  Coefficient  (8) 

with  permeability  (K)  and  porosity  (0) 
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Figure  13 
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Figure  14 
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A  plot  of  log  f  vs  log  F  was  also  attempted  but 

o 

it  showed  a  much  wider  scatter  than  the  plot  of  log  fs  vs 
log  0.  The  latter  was  therefore  chosen  in  developing  the 
proposed  correlations „ 

Several  other  plots  were  drawn  which  indicated  the 
interrelationships  among  ft,  0,  K,  F  and  t.  Some  of  these  are 
shown  in  Figures  7,  8,  9,  10a  and  10b » 

The  formation  resistivity  factor  data  were  obtained 
using  both  the  two-electrode  and  the  four-electrode  methods. 
Dunlap (2 5)  reported  that  the  contact  resistances  may  develop 
at  the  current  electrodes  due  to  either  poor  connections  betw<een 
the  core  end  areas  and  the  electrodes,  or  low  core  saturations. 
The  development  of  the  contact  resistances  was  checked  by 
employing  the  check  proposed  by  Dunlap (25) .  The  electrical 
resistance  data  obtained  on  all  twenty  cores  is  summarized 
in  Table  3.  Data  on  the  contact  check  resistance  shows  that 
there  was  negligible  contact  resistance  and  therefore  the 
formation  resistivity  factor  values  are  reliable. 

Several  materials  were  tried  to  make  contact 
between  the  core  and  the  current  electrodes  so  as  to  allow 
passage  of  electric  current  through  the  core,  without  any 
appreciable  development  of  contact  resistance.  Brine 
saturated  kleenex,  sponge,  Chamois  leather,  felt  pad  and 
drilling  mud  were  tried  but  none  was  found  satisfactory. 

Finally,  a  felt  pad  threaded  with  thin  silver  wire  to  form  a 


i 


aoi  o.io 

d> ,  (fraction) 


46 


Figure  10a  Figure  10b 


4 


3 


2 


1 


0 


-1 


-2 


-3 


-4 
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% 

* 

o 

o 


(mds) 


Table  3 


Formation 

Resistivity  Factor  and 

Tortuosity  Dat 

a 

T 

**  * 

Core 
No . 

Pb 

ohms-cm 

Rb 

ohms** *** 

Rc4 

ohms 

Rc2 

ohms 

Rcc 

ohms 

F 

0 

% 

3 

13.53 

10.15 

80.3 

300 

150 

29.56 

18.82 

2.36 

4 

13.53 

9.96 

62.3 

181.5 

90.5 

18.23 

21.11 

1.96 

5 

13.53 

10.00 

78.8 

230 

117 

23.0 

17.16 

1.99 

7 

13.53 

10.23 

81.3 

243.5 

123 

23.81 

13.60 

1.80 

8 

13.53 

10.09 

47.2 

17  5 

86 

17.34 

19.39 

1.83 

10 

13.53 

9.50 

* 

181 

90 

19.05 

19.96 

1.95 

11 

13.53 

8.82 

* 

145 

73 

16.44 

20.39 

1.83 

12 

13.53 

9.30 

* 

300 

145 

32.30 

9.82 

1.78 

13 

13.53 

9.65 

* 

790 

392 

82.0 

10.28 

2.90 

14 

13.53 

9.73 

* 

283 

143 

29.07 

18.85 

2.34 

4X 

13.53 

10.50 

45.0 

135 

65 

12.79 

22.18 

1.68 

20 

14.1 

10.73 

* 

1200 

587.5 

111.83 

2.30 

1.60 

22 

13.53 

10.48 

49.0 

142 

72 

13.55 

24.33 

1.82 

24 

13.53 

10.69 

45.9 

170 

84.5 

15.91 

20.91 

1.82 

31 

14.1 

10.93 

160 

530 

262 

48.47 

5.84 

1.68 

38 

13.53 

10.67 

37.9 

136 

68.0 

12.72 

23.36 

1.72 

45 

13.53 

10.55 

39.0 

157.5 

78.2 

14.92 

21.55 

1.79 

46 

13.53 

10.58 

50.8 

171 

82.5 

16.16 

18.10 

1.71 

53 

13.53 

10.69 

40.5 

145 

72 

13.57 

22.06 

1.73 

61 

13.53 

10.81 

33.8 

127 

64 

11.75 

27.43 

1.80 

*  Core  dimensions  such  that  the  probes  could  not  make  contact. 

**  Resistance  of  the  same  ( congruent ) volume  of  brine  solution. 

***  (t)^  -  F0  =  Le/L  Reference  (5) 

=  resistivity  of  brine 
R-^  =  resistance  of  brine 

Rc^  -  resistance  measurement  by  4  electrode  method. 

Rc2  ”  resistance  measurement  by  2  electrode  method. 

Rcc  =  resistance  of  core  contact  resistance  check. 

F  =  formation  resistivity  factor 
0  =  porosity 


t  =  tortuosity 
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continuous  mesh  on  both  sides,  and  saturated  with  brine, was 
found  to  be  the  most  suitable  of  all  the  techniques  proposed 
so  far. 

Cornell (14)  obtained  data  on  formation  resistivity 
factor  by  employing  an  apparatus  in  which  the  core  sample  was 
placed  vertically  between  two  platinum  electrodes.  During  this 
investigation  it  was  found  that  the  vertical  position 
increased  the  drainage  of  brine  from  the  core  and  thus  affected 
the  resistance  measurements.  Also  Cornell! 14)  did  not  in¬ 
corporate  any  check  for  detecting  the  presence  of  contact 
check  resistance.  These  two  points  may  account  for  the 
scatter  observed  in  Figures  7  and  10a  showing  the  data 
obtained  by  Cornell(14),  Hamilton(7)  and  the  author. 

There  is  some  disagreement  over  the  quantitative 
relationship  between  the  tortuosity  and  the  formation 
resistivity  factor.  Cornell(14),  Street(27)  and  Wyllie  and 
Spangler (17)  have  suggested  different  models  thus  resulting 
in  different  expressions. 

In  order  to  compare  the  tortuosity  data,  obtained 
in  this  study  with  Cornell's  data,  a  model  for  tortuosity 
similar  to  that  suggested  by  Cornell  and  Katz(15)  was  used. 

Several  plots  were  tried  to  correlate  the  permeability 
with  different  powers  of  tortuosity.  A  plot  of  log  (K)  *  versus 
log  r  was  found  to  show  the  least  scatter  of  data. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

1.  A  general  correlation  has  been  developed 
relating  |3  to  porosity,  permeability  and  the  formation  resis¬ 
tivity  factor.  Unlike  other  correlations,  this  is  dimensionally 
consistent  and  is  not  based  on  the  Kozeny  model  of  porous  media. 

2 .  The  experimental  data  indicate  that  the  reservoir 
parameters,  i.e.  0,  K,  F  and  t  are  interrelated. 

3.  The  data  on  Klinkenberg  permeability  were  checked 
with  the  permeability  values  obtained  by  evaluating  a  values. 
Good  agreement  v/as  observed. 

4.  The  resistance  measurements  by  the  two  electrode 
and  the  four  electrode  methods  are  reliable  only  if  no  contact 
resistance  is  allowed  to  develop  at  the  current  electrode. 

5.  A  brine  saturated  felt  pad,  threaded  with  thin 
silver  wire  so  as  to  form  a  continuous  mesh  on  both  sides,  was 
found  more  suitable  as  a  contact  material  between  the  core  and 
the  electrode  than  any  other  method  reported  in  the  literature. 

6.  Anomalous  Klinkenberg  permeability  data  would  be 
obtained  if  proper  care  is  not  exercized  to  eliminate  the 
effect  of  slippage. 

Although  an  improved  correlation  between  the  inertial 
resistance  coefficient  (3  and  the  reservoir  rock  properties  is 
presented  herein,  it  should  not  be  considered  as  a  terminal 
attempt .  The  following  recommendations  are  therefore  made  for 
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future  work: 

1.  The  investigation  of  other  measurable  properties 

of  porous  media  for  their  influence  on  (3,  K  and  0  is  recommended. 
Some  of  those  properties  are  the  internal  surface  area,  tor¬ 
tuosity  by  a  means  independent  of  formation  resistivity  factor, 
surface  roughness,  pore  size  distribution,  etc. 

2.  For  obtaining  data  on  permeability  of  very  tight 
rock  samples,  it  is  recommended  to  further  investigate  the 
prospects  of  using  Selas  porcelain  microfilter  discs  and 
candles  as  a  volume  measuring  device. 

3.  The  effectiveness  of  applying  back  pressure  to 
restrict  the  flow  to  the  viscous  region  should  further  be 
studied.  The  Klinkenberg  permeability  data  should  be  obtained 
on  several  cores,  at  constant  mean  pressures  but  with  different 
flow  rates,  in  order  to  determine  those  flow  rates  that  will 
have  to  be  achieved,  before  the  assumption  that  the  data 
points  fit  the  straight  line(ll)  is  necessarily  true. 

4.  The  problem  of  finding  a  suitable  contact  material 
to  carry  the  electrical  current  between  the  current  electrodes 
and  the  core  should  further  be  examined.  To  minimize  errors 
the  resistance  measurement  cell  should  be  enclosed  in  a  con¬ 
stant  temperature  and  humidity  cabinet. 
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NOMENCLATURE 


A 


D 


F 


f 


s 


G,  Gs 
Uc 


K 

K1 


Q,  q 


Cross  sectional  area  of  porous  sample,  square  inches 
Diameter  of  porous  sample,  inches 
Effective  pore  diameter  of  pore  structure,  (feet) 
Electrical  resistivity  factor 

A  factor  defined  by  grouping  and  K  to  develop 
the  proposed  general  correlation 

Mass  velocity  (pound  mass/square  foot  second) 

Conversion  factor  (32.17  pounds  mass  x  foot/ 
pounds  force  x  second  squared) 

Permeability  of  a  porous  medium  to  a  fluid, 
millidarcies  or  (ft^) 

Dimensionless  geometrical  constant  in  the  "laminar" 
term  of  flow  equation  -  reference  (14) 

Dimensionless  geometrical  constant  in  the 
"turbulence"  term  of  flow  equation  developed  by 
Cornell  -  reference  (14) 

Apparent  permeability,  millidarcies 

2 

Klinkenberg  permeability  (ft  ) 

Length  of  porous  sample,  inches 

Fictitious  average  path  length  through  a  porous 
medium  (ft) 

Molecular  weight  of  a  gas 

Upstream  Pressure,  (pounds  force/sq. inch)  absolute 
Downstream  pressure,  (pounds  force/sq . inch)  absolute 

Mean  pressure,  (pounds  force/sq.  inch)  absolute 
Flow  rate,  cubic  feet/hour 
Flow  rate  @  standard  conditions 


52 


R 


T 


v 

x 

Z 


a 

P 

A 

P* 

P 

pb 


Gas  constant  (1543  pounds  force  x  cubic  feet/ 
square  feet  x  pound  mole  x  °R) 

Resistance  of  core, ohm 

.  2 

Internal  surface  area  of  a  porous  medium,  ft 
Absolute  temperature,  °R 
Velocity,  (feet/second) 

Length  of  porous  medium 
Compressibility  factor  of  a  gas 
Viscous  resistance  coefficient,  (ft  ) 
Inertial  resistance  coefficient,  (ft~^) 
Difference 

Viscosity  of  gas  (pounds  mass/foot  second) 
Density  of  gas,  (pounds  mass/cubic  feet) 
Resistivity  of  brine,  ohm  -cm 
Tortuosity 
Porosity 


0 
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Brine  Saturation  Data 


57 


G 

0 

•H 

-P 

rd 

p 

G 

•P 

rd 

w 


oirio  <yw.-icor-'r^cMcoiDOcocoincr»cOrH 


id  vd  in 

(Ti  CP\  (T>  (T> 


r^covo^  'd*vor-<Nr^r^cMr^r'MDr^co 

C?i(T'0^0^CTiC?i0^C?>C^C^0^C?'0^G^Cr>0^ 


0 

u 

o 

PM 


o 


£ 


G  O 

r— I  (J 

o 


> 


CO  LO 
<T>  ^ 
LO  <T> 


cm  in 
o  o 

CM  ID 


id  m 

"sf  CM 
ID  CO 


CO  CM 
CM 

ID  CO 


CM  rH 
O  O 

O  r» 


r-(MOcocooofoincri 

rH^flDCMOinONr^OlD 


cocncocMcocorOfHCMcn^O'd'^iH'^cocDinin 


4H 

0 

0 

£ 

G 

rH 

o 

> 


TS 
0 
0  -P 
G  -P  U 
•H  -H  U 


0 


OrH^inC0^r^lDlD(T>fnC0OCMCMrHf0^cr>CM 

'sDcor-r^^tv£»co^<^r^incocN>inooOiHroco 

^0O^mrv^r^0inon^i-io^00(^in 


mnncNcnncoHr-m^o^^H^nn^in 


>i 

-P 

•H 

0 

G 

0 

Q 


0  U 

CTi 

<T> 

cr> 

C?1 

CT>  <T> 

C7> 

O'! 

O'* 

CT> 

cr> 

<T> 

cr> 

<T> 

<y> 

<T> 

<T> 

<r> 

o>  cr> 

GO 

CTi 

CP> 

c?> 

CTi 

cri  cr\ 

o> 

O'* 

<T> 

cr> 

<T> 

0^ 

cr> 

<g 

CT> 

<T> 

(T»  CT> 

4-1 

•H  \ 

CM 

CM 

CM 

Csl 

CM  CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM  CM 

0 

*1 

O 

O 

o 

O 

O  O 

O 

O 

O 

O 

O 

O 

O 

O 

o 

o 

O 

O 

O  O 

rH 

fH 

rH  f — 1 

i — 1  rH 

rH 

rH 

rH 

«H 

rH 

rH 

rH 

rH 

i — 1  rH 

rH 

rH 

rH  rH 

4H  Td 
0  0  ^ 
0  -P  0 
-P  G  -P  £ 
,G  -H  -H  U> 
t^P  £  — 
■H  ,0  rO 
0  0 
£ 


min<yiOrHincr)cMr'OCMiDcnrHOiDr^rHr>-<?i 

COinin<T>Td,lDHC0CMlD'd,Cr>'d,lDCM<MlDCMlDC0 

iD,d,iDTd,inr'Cr»cr>incor'Cr>cMr'CorOr-icMCO,:t 

in<T»rHiniDcoinr^<j'iDiHcniDcMOino^r^or^ 


ncnrocMnncOr-iHo^o^^H^roninin 


4M 

"0 

0 

0 

CO 

rH 

(7> 

i H 

cr> 

CM 

r- 

-p 

0 

rH 

CO 

CM 

r' 

r- 

-P 

G 

P 

0 

1 — 1 

CM 

CG 

in 

CO 

O 

in 

rG 

P 

0 

£ 

o 

1 - 1 

o 

r- 

LO 

r- 

tn 

G 

O 

& 

•H 

-P 

** — ^ 

in 

CO 

r- 

CO 

o 

0 

rd 

& 

0 

o 

o 

CM 

ID 

r- 

o 

CO 

i — 1 

cr> 

r- 

in 

rH 

rH 

cr> 

CM 

(G 

CO 

ID 

o 

CO 

CM 

CT> 

rH 

in 

CO 

(T> 

in 

CO 

rH 

r- 

CT> 

CT> 

o 

CO 

CM 

CO 

ID 

in 

in 

ID 

CO 

ID 

CM 

CO 

rH 

in 

r- 

CO 

1 - 1 

in 

■'t 

CO 

•sf 

'si1 

ip 

0 

0 

CO 

ID 

o 

rH 

CO 

r-* 

00 

r- 

O 

O 

CO 

CO 

CO 

r- 

CO 

r- 

CTi 

ID 

in 

P 

CO 

in 

CO 

CO 

CO 

r~ 

ID 

CM 

ID 

rH 

rH 

CM 

o 

CO 

c?> 

ID 

in 

ID 

rH 

in 

-P 

o^. 

in 

r- 

CM 

o 

CO 

CM 

ID 

ID 

in 

ID 

ID 

CM 

CO 

1 — 1 

in 

in 

rH 

o 

,G 

0  0 

rH 

<r> 

CM 

rH 

ID 

rH 

O 

rH 

O 

CM 

in 

CO 

p" 

CO 

f'- 

<T> 

cr> 

in 

00 

tr> 

£ 

•H 

>i  &> 

rH 

cr> 

o 

m 

o 

o 

t"- 

in 

ID 

CO 

CO 

00 

r- 

o 

ID 

CO 

r- 

rH 

<r> 

r- 

0 

P  ^ 

CO 

CO 

'vf 

CO 

CO 

CO 

CO 

CO 

CO 

& 

0Z5 

0 

p  • 

O  O 
U  £ 


CO^inr-'COOiHCMC’O^XOCM^trHCOinvDCOrH 

HHHHH^fN)<MCMrm^^in0 


V  • 

1 

2 

3 

4 

5 

6 

7 

G 

9 

10 

11 

12 

13 


58 


Table  5 

Viscous  Flow  Data 
Core  Number  4X 


Pressure,  psia 

1 


P1 

P2 

Pm 

(psia) 

(Pi2  -  p22) 
(psia) 2 

Tf 

°R 

Qc 

SCFH 

K 

(md) 

14.97 

13.465 

0.0703 

42.88 

537.0 

0.058 

65  #2 

17.41 

13.479 

0.0648 

121.29 

536.2 

0.158 

62.2 

19.90 

13.612 

0.0597 

210.56 

534.0 

0.275 

62.2 

22.83 

13.516 

0.0550 

338.39 

535.8 

0.434 

61.4 

26.42 

13.550 

0.0500 

514.52 

535.7 

0.669 

62.2 

30.65 

13.597 

0.0452 

754.42 

535.7 

0.967 

61.4 

36.52 

13.684 

0.0398 

1146.24 

535.6 

1.446 

60.4 

42.15 

14.676 

0.0352 

1561.24 

535.7 

1.957 

60.0 

47  *  45 

19.159 

0.0300 

1884.2 

531.2 

2.39 

60.3 

53.47 

26.685 

0.0250 

2146.8 

531.3 

2  .689 

59.5 

64.55 

35.71 

0.0200 

2891.4 

531.6 

3.551 

58.4 

80.29 

52  .60 

0.0151 

3679.1 

532  .9 

4.409 

57.2 

118.4 

92.21 

0.0095 

5525.1 

533.0 

6.380 

55.2 
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Table  7 

Viscous  Flow  Data 
Core  Number  20 


Pressure,  psia 

P1  P2 

1 

Pm 

(psia)"1 

„  2  _  2 
P1  "  P2 
(psia) 2 

Tf 

°R 

Qc 

SCFH 

260.3 

13.307 

0.0073 

67577 

532.2 

0.000399 

316.2 

13.298 

0.0061 

99804 

532.2 

0.000532 

37  5.1 

13.298 

0.0052 

140522 

532.2 

0.000715 

437.3 

13.301 

0.0044 

191054 

532  .2 

0.000888 

488.8 

13.300 

0.004 

238749 

532.6 

0.00105 

586.8 

13.300 

0.0033 

344157 

532  .7 

0.00142 

693.3 

13.300 

0.0028 

480488 

532  .9 

0.00187 

839.3 

13.302 

0.0024 

704251 

533.1 

0.00250 

985.7 

13.304 

0.0020 

971435 

533.3 

0.00328 

1122.3 

13.309 

0.0018 

1259394 

533.7 

0.00403 

1391.9 

13.311 

0.0014 

1937231 

534.0 

0.00581 

1543.1 

13.309 

0.0013 

2381005 

534.1 

0.00677 

1659.3 

13.307 

0.0012 

2753119 

534.3 

0.00779 

K 

(md) 

0.000275 

0.000249 

0.000238 

0.000218 

0.000205 

0.000194 

0.000184 

0.000168 

0.000161 

0.000153 

0.000144 

0.000138 

0.000135 
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Table  9 


Viscous  Flow  Data 
Core  Number  22 


Pressure,  psia 


p 

1 

P 

2 

1 

pm 

(psia) 

2  2 
v  -  y 

(psia)2 

Tf 

°R 

Qc 

SCFH 

K 

(md) 

16.45 

13.809 

0.0661 

79.81 

532.3 

1.716 

1020 

18.38 

14.582 

0.0607 

125.30 

532.3 

2  .673 

1012 

19.29 

15.393 

0.0577 

135.08 

532.7 

2  .749 

956.5 

19.83 

16.342 

0.0553 

126.05 

532  .8 

2.593 

977.2 

21.68 

18.214 

0.0501 

138.45 

533.0 

2.807 

964 

23.87 

20.314 

0.0453 

157 .21 

533.1 

3.134 

948 

26.59 

23.608 

0.0398 

149.58 

533.3 

2.982 

940 

30.58 

27.274 

0.0346 

191.45 

533.5 

3.696 

918.5 

34.82 

31.93 

0.0300 

193.05 

533.7 

3.697 

912 

41.63 

38.73 

0.0249 

233.53 

533.8 

4.329 

883 

50.87 

48.50 

0.0201 

236.0 

534.0 

4.305 

87  0 

67.54 

65.42 

0.0150 

281.87 

534.0 

4.683 

793 

101.27 

98.42 

0.0100 

569.09 

537  .9 

8.65 

7  32 

101.25 

98.70 

0.0100 

509.87 

534.0 

7.31 

692 

Visco- Inertial  Flow  Data 
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Table  11 

Viscous  Flow  Rate 
Core  Number  24 

Run  Pressure,  psia 


• . 

pi 

P2 

JL 

m 

(psia)"1 

2  2 
P1  -  P2 

(psia)2 

Tf 

°R 

Qc 

SCFH 

K 

(md) 

1 

19.66 

13.660 

0.0600 

199.96 

537 .3 

2.067 

505 

2 

21.50 

14.777 

0.0551 

243.72 

533.0 

2.477 

492  .7 

3 

23.43 

16.398 

0.0502 

279.88 

533.3 

2.772 

480.7 

4 

25.54 

18.464 

0.0455 

311.47 

533.5 

3.011 

469.3 

5 

28.55 

21.405 

0.0400 

356.76 

533.8 

3.360 

457.6 

6 

31.64 

25.481 

0.0350 

351.94 

534.1 

3.308 

456.9 

7 

36.60 

29.995 

0.0300 

439.71 

534.3 

3.965 

438.5 

8 

52.60 

47  .27 

0.0200 

532 .67 

534.9 

4.500 

412.1 

9 

101.25 

98.35 

0.0100 

578.84 

535.6 

4.601 

388.5 
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Table  13 

Viscous  Flow  Data 
Core  Number  31 


Run  Pressure,  psia 


No. 

pi 

P2 

1 

P™ 

m 

(psia) 

2  2 

P  -  P 

1  *  2 

(psia) ^ 

\ 

Tf 

°R 

Qc 

SCFD 

K 

(md) 

1 

83.5 

13.441 

0.0206 

6972 

533.8 

0.00121 

0.00826 

2 

103.4 

13.441 

0.0171 

10519 

534.2 

0.00164 

0.00744 

3 

116.5 

13.483 

0.0154 

13387 

534.9 

0.00195 

0.00697 

4 

129.4 

13.441 

0.0140 

16574 

534.4 

0.00234 

0.00674 

5 

163.6 

13.489 

0.0113 

26570 

535 

0.00353 

0.00636 

6 

199.5 

13.504 

0.0094 

39620 

535.1 

0.00476 

0.00576 

7 

223.5 

13.528 

0.0084 

49782 

531.0 

0.00588 

0.00562 

8 

273.5 

13.533 

0.0070 

74634 

531.5 

0.00832 

0.00532 

9 

317.6 

13.550 

0.0059 

100708 

531.9 

0.0109 

0.00517 

10 

419.6 

13.554 

0.0046 

17  5859 

532.5 

0.0178 

0.00483 

11 

535.5 

13.554 

0.0036 

286619 

533.0 

0.0267 

0.00448 

12 

533.5 

13.545 

0.0037 

284471 

534.0 

0.0261 

0.00443 

13 

603.5 

13.547 

0.0032 

364064 

534.3 

0.0323 

0.00428 

14 

703.5 

13.492 

0.0028 

494688 

537.9 

0.0409 

0.00403 

15 

808.5 

13.497 

0.0024 

653438 

537  .9 

0.0524 

0.00392 

16 

929.5 

13.504 

0.0021 

863821 

538.0 

0.0664 

0.00377 

17 

1065.5 

13.507 

0.0019 

1135133 

538.0 

0.843 

0.00366 

18 

1213.5 

13.495 

0.0016 

1472364 

533.7 

0.111 

0.00369 

19 

1382.4 

13.511 

0.0014 

1910944 

533.7 

0.139 

0.00357 

20 

1573.4 

13.531 

0.0013 

2475521 

533.9 

0.171 

0.00344 

' 

Viscous- Inertial  Flow  Data 
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Table  15 

Viscous  Flow  Data 
Coro  Number  J38 _ 


Run  P  r  t !  s  s  u  r  o ,  p  s  .1  a 


o . 

pi 

P2 

1 

pm 

( p  s  i  a )  ~ 

2  2 
P1  -  P2 

(psia) 2 

Tf 

°R 

0C 

GCFII 

K 

(md) 

1 

16.94 

14.072 

0.0645 

89.1 

536.9 

2.272 

1244 

2 

19.06 

14.700 

0.0592 

147  .2 

536.8 

3.717 

1230 

3 

19.32 

15.308 

0.0578 

138.9 

536.5 

3.486 

1220 

4 

20.51 

16.316 

0.0543 

154.5 

536.5 

3.840 

1210 

5 

22.01 

18.015 

0.0500 

159.9 

536.5 

3.801 

1190 

6 

24.01 

21.124 

0.0443 

130.1 

536.5 

3.109 

1165 

7 

26.13 

23.74 

0.040 

119.0 

536.5 

2  .879 

1178 

8 

29.96 

27.45 

0.035 

143.8 

533.3 

3.411 

1149 

9 

34.34 

32.50 

0.030 

122  .8 

533.5 

2.943 

1161 

10 

41.06 

39.03 

0.02  5 

162.5 

533.7 

3.706 

1106 

11 

135.45 

134.36 

0.0074 

294.0 

535.0 

6.174 

1027 
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Table  17 

Viscous  Flow  Data. 

Core  Number  45 


Pressure,  psia 


P1 

P2 

1 

pm 

(psia)  1 

2  2 

P  -  P 

F1  *2 

(psia)2 

Tf 

°R 

SCFH 

K 

(md) 

14.77 

13.438 

0.0709 

37.7 

535.1 

0.518 

661.5 

16.54 

13.570 

0.0664 

89.4 

535.1 

1.118 

601.1 

16.48 

13.979 

0.0657 

76.1 

535.2 

1.010 

638.9 

18.52 

14.709 

0.0602 

126.5 

535.2 

1.608 

611.7 

20.25 

15.407 

0.0561 

172.6 

535.7 

2.138 

596.5 

22.44 

17.594 

0.0500 

194.1 

535.7 

2.350 

582.9 

23.64 

20.768 

0.0499 

127.5 

535.7 

1.629 

570.0 

25.60 

23.810 

0.0403 

88.6 

535.8 

1.261 

558.0 

29.15 

27 .476 

0.0352 

94.5 

535.8 

1.384 

545.0 

33.91 

32.505 

0.0300 

93.5 

536.0 

1.354 

535.1 

40.74 

39.448 

0.0249 

103.8 

536.0 

1.525 

518.0 

52.11 

47.252 

0.0201 

482  .9 

537 .5 

5.023 

505.2 

101.24 

98.292 

0.0100 

588.6 

538.3 

5.794 

479.5 
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Table  19 

Viscous  Flow  Data 
Core  Number  46 


Pressure,  psia 


P1 

P2 

_1_ 

pm 

(psia)"1 

2  2 
P1  “  P2 
(psia) 2 

Tf 

°R 

Qc 

SCFH 

K 

(m<3) 

15.09 

13.53 

0.070 

44.7 

540.1 

0.085 

93.3 

16.88 

13.53 

0.066 

105.1 

540.0 

0.195 

94.0 

19.32 

14.06 

0.060 

175.4 

540.0 

0.333 

93.0 

19.88 

14.51 

0.058 

184.6 

531.3 

0.354 

92.5 

21.41 

14.84 

0.055 

238.2 

532.0 

0.455 

92.0 

22.37 

18.15 

0.049 

170.8 

532.2 

0.322 

91.0 

24.77 

20.08 

0.045 

210.4 

532.3 

0.390 

89.5 

27.10 

22.63 

0.040 

222.3 

532.3 

0.411 

89.3 

30.20 

26.38 

0.035 

216.3 

532.7 

0.395 

88.2 

34.63 

31.89 

0.030 

182.2 

533.0 

0.329 

87.4 

41.44 

38.76 

0.02  5 

214.6 

533.3 

0.386 

87.0 

51.49 

48.31 

0.020 

317.4 

533.7 

0.551 

84.1 

100.8 

98.6 

0.010 

446.5 

533.8 

0 . 660 

71.7 
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Table  21 

Viscous  Flow  Data* 
Core  Number  53 


Run  Pressure,  psia 


1  • 

pi 

P2 

1 

pm 

(psia)”1 

2  2 
P1  -  P2 

(psia)2 

Tf 

°R 

Qc 

SCFH 

K 

(md) 

1 

14.56 

13.509 

0.0713 

29.50 

537.0 

0.282 

469.0 

2 

17  .90 

13.639 

0.0634 

134.50 

537  .8 

1.263 

460.9 

3 

21.50 

13.884 

0.0565 

269.66 

538.0 

2.472 

450.0 

4 

23.77 

18.248 

0.0476 

231.88 

537 .0 

2  .082 

440.6 

5 

29.03 

22.07 

0.0391 

355.42 

538.3 

3.090 

425.1 

6 

37.63 

30.60 

0.0293 

479.66 

536.1 

4.037 

411.0 

7 

52.25 

46.45 

0.0203 

572.48 

536.0 

4.708 

402.1 

8 

101.45 

98.3 

0.0100 

623.30 

536.5 

4.814 

37  5.0 

9 

136.50 

133.5 

0.0074 

810.00 

536.7 

5.985 

360.0 

*  The  data  were  obtained  by  applying  back  pressure. 
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Table  22 

Anomalous  Data  for  Klinkenberg  Permeability* 

Core  Number  53 


Run  Pressure,  psia 

No.  1 


P 

1 

P 

2 

P 

m 

(psia)"1 

2  2 

P1  -  P2 

(psia)2 

Tf 

°R 

Qc 

SCFH 

K 

(md) 

1 

15.36 

13.314 

0.0698 

59.0 

535.5 

0.54 

447 

2 

18.87 

13.402 

0.0620 

176.0 

535.5 

1.59 

439.5 

3 

22.33 

13.544 

0.0558 

314.1 

535.3 

2  .71 

418.5 

4 

26.14 

13.670 

0.0503 

498.0 

533.5 

4.09 

397.3 

5 

30.13 

14.025 

0.0453 

712  .0 

535.3 

5.50 

374.0 

6 

34.37 

14.34 

0.0411 

974.0 

533.8 

7.10 

352.0 

7 

51.25 

16.828 

0.0294 

2351 

535.3 

13.90 

288.0 

8 

88.05 

23.79 

0.0197 

7213 

538.3 

30.60 

206.0 

9 

114.8 

31.33 

0.0136 

12216 

535.2 

43.50 

172.5 

10 

140.3 

39.34 

0.0112 

18450 

534.8 

56.16 

148.0 

11 

167.3 

48.2  5 

0.0093 

25670 

534.3 

64.8 

122.6 

*  Calculated  without  back  pressure. 
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i — I  rH  rH  i — I 


<N  oo  co  m  in  id  incN^c^cNincor^oocri 
ocoa^a>incoinr-incriOiDcniDinc40 
ocor'Or^n'inOiHcocNinrHroococo 

^rOr-lDrHCNOiD^lD^rHr^CNCOOO 

rH(NrHCOrH'2’rHlDCNCr>CNCOCnr- 


coinoOrHrHcor^roocnfncnr^cNOco 

r-r-coooco^LncomiDCOiDcocr)inr--^t 

cocncncncncocococncocorncocncocnrn 

ininininininminininininininininin 


corHiH^in  in  (n  co 

iDiHlDOOl^^rHCOOCOCMCN  CO  OlD 
CNCN^fOCOOlDOO^rHCNCOlDOiniDrH 

rHiHCMinuo^c^incoinoiDOr'CncniD 

rH  CN  cn  ^  in 


^coc^cNOOcncocoin^tiDc^iDcocri^ 

iDincor'iniDcNCOCO^iDcNp-cncoocn 

cocoro^'inomLniDiDCOcocor-rHCOc^ 

r!HHHHfOHfOH^H0(N(TifOf\lfO 


rH^rHo^cococor^Lnin  nmomH^ 
o^cr>mcnr-iDOOiDCNCN^oo^tr^^tco 

r^r^(HCOfnr^OfOrHcNrHOCOCN^tcn<7> 

rHrHCNCNcooo^'itLniniDn'COOrHmco 

rH  i — I  rH  rH 


rHiNco^mior^coo^OrHcNC'Q'itinior--- 

i — I  rH  i — I  rH  i — I  H  rH  rH 
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Table  24 


Viscous  Flow  Data 
Core  Number  61 


Run  Pressure,  psia 

No.  1 


P1 

P2 

Pm 

m 

(psia) 

2  2 

P  -  P 

F1  F2 

(psia)2 

Tf 

°R 

SCFH 

K 

(md) 

1 

19.24 

14.46 

0.0594 

161.2 

538.2 

4.545 

1395 

2 

21.04 

15.01 

0.0555 

217.1 

538.2 

5.896 

1344 

3 

23.16 

16.89 

0.0999 

251.3 

538.2 

6.656 

1311 

4 

2  5.37 

19.40 

0.0447 

267.5 

538.2 

7.010 

1298 

5 

27.60 

22.38 

0.0400 

260.9 

538.2 

6.837 

1297 

6 

31.50 

25.70 

0.0350 

331.9 

538.2 

8.349 

1245 

7 

35.93 

31.03 

0.0299 

327.8 

538.2 

8.258 

1248 

8 

41.84 

37  .92 

0.0251 

313.2 

538.2 

7.906 

1250 

9 

51.98 

47  .86 

0.0200 

511.3 

538.2 

10.177 

1226 

10 

100.85 

98.81 

0.0100 

409.3 

538.2 

9.284 

1126 
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3°-w(sd) 


Figure  2-A.  Flow  Data  Plot  to  Calculate  a  and  3 
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Figure  3-A.  Klirvkenberg  Permeability  Plot 


Fig.  4— A.  Flow  Data  Plot  to  Calculate  a  and  p 
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Figure  6-A.  Flow  Data  Plot  to  Calculate  a  and  (3 
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Figure  7-A.  Klinkenberg  Permeability  Plot 


Figure  8-A.  Flow  Data  Plot  to  Calculate  a  and  0 


83 


Figure  9-A.  Klinkenberg  Permeability  Plot 


Figure  10-A.  Flow  Data  Plot  to  Calculate  a  and  p 
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Figure  12- A.  Flow  Data  Plot  to  Calculate  a  and  £ 
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Figure  14-A.  Flow  Data  Plot  to  nslmiaf 
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Figure  16— X.  Flow  Data  Plot  to  Calculate 
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Figure  18-A.  Flow  Data  Plot  to  Calculate 
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EQUATIONS 


2  Qc  Pb 

(pi  +  P2) 


Qm  =  flow  rate  at  mean 

pressure  and  flowing 
temperature 

Qc  =  flow  rate  at  standard 
conditions,  14.65  psia 
and  60 °F 

Pb  =  pressure  base  =  14.65  psia 

=  average  flowing  tem¬ 
perature  in  °R 

Tb  =  temperature  base  =  60  °F 

P-^  =  upstream  pressure 

P 2  =  downstream  pressure 


Qm 

A  pi  -  p2 


Mb  ^  Tf 
A  Tb 


K  =  darcies 
a 

[i  =  CP 

Q  --  cc/sec 

L  =  cm 

P  =  atm 

A  =  cm2 

°R 


T 


* 
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K. 


=  0.25573 


Qc  Tf  \x  h 
(AP2)  A 


K. 


=  millidarcies 


_  o 


R 


(1 

L 

A 

P 

CL 


=  micropoise 


=  m 


,  2 
in 


psia 

SCFH 


M  g 


2  (ZRT)  \x  L 


=  52.0956  x  106 


lbm 


ZT^pJj  lb^  sec 


j  ”  G 


Qcp 


s 


04 


Qcp  lb, 


m 


A  ft  sec 


A(P2)  A(P2 )  lbf2  sec 

- —  =  207  36  — — - 


G  lbm  ft' 


Z 

Tf 

|l 

L 

Q 

P 

A 

P 


=  compressibility 


=  °R 


=  micropoise 


=  m . 


=  SCFH 

=  lbm/ft3 
=  in2 


psia 


G  5  G  _i 

-  =  148.8  x  10D  -  ft  ± 

\i  p. 


M- 


lb. 


ISU 


ft  sec 
micropoise 


meter  correction  factor  x 


x  35.495 


+  P 

meter  m3 


Q  =  CFH 

=  meter  inlet  pressure 

in  psia 

Tm  =  meter  inlet  temperature 


.et  (psig) 


partial  pressure  of  water 
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SAMPLE  CALCULATIONS 

Klinkenberg  Permeability  Data 

Core  No.  53  -  Taber  sandstone 

Area  =  0.7775  in^  Length  =  1.5593  in. 

Posority  =  22.06% 

Q  =  flow  rate  =  2.7975  CFH 

Temperature  =  535. 5 °R 

Pressure  atmospheric  =  69-87  cm  Hg  =  13.464 

Flowing  fluid  -  nitrogen 

Average  flowing  temperature  *  538.0  °R 

fl  (avg.  temp,  and  pressure)  =  177.875  micropoise 


P1 

— 

21 . 504  psia 

P2 

= 

13.884  psia 

pmeter 

= 

13.464  psia 

p 

base 

— 

14.65  psia 

m 

base 

520  °R 

A(P2) 

— 

269.66  (psia) 

water  differential  on  the  manometer  =  0 


Calculations 


1.  Q 


SCFH 


, .  _  p meter  Tbase 

0/  .  v  x  meter  correction  factor  x  — -  x  - 

-  (meter )  *P  p 


'meter 


base 


QSCFH 


2.4715 


. 
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2 .  K. 


0.25573 


QC  Tf  L 

A  (P2)  A 


-  450  mds 


3.  =  0.04 


Qrp 


9.362  x  10" 3  lbm/ft2 


Q 

4.  —  =  148.8  x  105  x  — 


G< 

\l 


=  783.17  ft 


-1 


5. 


A(P2 ) 


=  20736 


A(P2) 


=  5.9726  x  108  lb  2-sec/lb  -ft2 

f  m 


6. 


M  gc 


2  (ZRT)  lib 


52.0956  x  10' 


Tf  |jl  Z  L 


3.4907  x  102  lb  /ft2-sec 

m 


M  g_  A(P2)  X1  _2 

7.  - i: - —  =  2.0849  X  10  (ft  ) 

2  ( ZRT)  faLG 


sec 


